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Abstract. Suppose that M is a compact Riemannian manifold with bound- 
ary and u is an L 2 -normalized Dirichlet eigenfunction with eigenvalue A. Let 
tji be its normal derivative at the boundary. Scaling considerations lead one 
to expect that the L 2 norm of ip will grow as A 1 / 2 as A — > oo. We prove an 
upper bound of the form < CA for any Riemannian manifold, and a lower 

bound cA < \\1pW2 provided that M has no trapped geodesies (see the main 
Theorem for a precise statement). Here c and C are positive constants that 
depend on M, but not on A. The proof of the upper bound is via a Rellich-type 
estimate and is rather simple, while the lower bound is proved via a positive 
commutator estimate. 



1. Introduction 

Let M be a bounded domain in W 1 with smooth boundary, or more generally a 
smooth compact Riemannian manifold with boundary dM = Y . Let H = —Am 
be minus the Dirichlet Laplacian on M. It is well known that H has discrete 
spectrum < Ai < A2 < A3 • • • — * 00. Let Uj be an L 2 -normalized eigenfunction 
corresponding to Xj , and let ipj be the normal derivative of Uj at the boundary. In 
this paper we consider the following question: do there exist constants c > and 
C < 00, depending on M but not on j, such that 

(1.1) cXj < Uj\\h(Y) <C\j ? 

Note that these inequalities exhibit the expected scaling; the square of the L 2 
norm of ipj involves two derivatives in space, which should scale as Xj in the high 
energy limit Xj — > 00. This question was posed by Ozawa in Q]. Using heat kernel 



techniques, he showed that an averaged version of (1.1) holds. More precisely, he 
showed that 

\(n/2) + l 

£ *s <»> - (4^r((„/2 )+ 2) + °< A '*" 2H1) - <°' ™ h y <= Y 

Xj < A 



This asymptotic formula (after integrating over Y) would be implied by (LI) in 
view of Weyl asymptotics for the Xj . 

In this paper, we prove the upper bound from (1.1) in general, and obtain a 



result for lower bounds. The main theorem is given later in the introduction. 
Before stating it, we give some simple examples which show in particular that the 
lower bound fails for some Riemannian manfolds with boundary. Additional, more 
elaborate examples are given in Section |^. 



A. H. is supported by an Australian Research Council Fellowship. T. T. is a Clay Prize Fellow 
and is supported by the Packard Foundation. 
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Example 1 — the disc. Let M = {x G M 2 | \x\ < a} for some a > 0. In this case 
we have an equality 

(i.2) / ^(9)de = ^i. 

J Si a 

This follows from Rellich's identity which also gives the same result in all 

dimensions. 

Example 2 — the rectangle. Let M — [0, a] x [0, b], where a < b. Then normalized 
eigenfunctions take the form 

/4 rmrx rimy , mils 2 ,nn-.2 
— sm sm — — , m,n integers, A = I I + ( — ) . 
ab a b a b 

The square of the L 2 norm of its normal derivative is 

1 mr rmrx 2 f b ran mry 2 \ 

(—sin ) dx + ( sin— — ) dy 

b a J a b J 

4 / / T17T \ 2 , TO7T N 2,\ 

From this we see that 

(i-3) \\<m\<-\ 

b a 

and by taking morn equal to one and sending the other integer to infinity, we see 
that these bounds are the best possible. Thus, for a square, we get equality as we 
do for a disc, but in general, the ratio between the best upper and lower bounds is 
the ratio of the side lengths. 

Example 3 — the cylinder. Let M = [0, a] xSj, the product of an interval with 
a circle of length b. Then the normalized eigenfunctions are 

— sm e'™ 9 ' 6 , m, n integers, A=( ) + ( — ) . 

ab a a b 

The square of the L 2 norm of its normal derivative is 

2 A ( b \™j™e/^ de 



ab 



a 

a a 



2 



In this case, the upper bound 

Mla<-A 
a 

holds, but no lower bound holds since we may hold m fixed and send n to infinity. 

Example 4 — the hemisphere. Let M be the hemisphere 

M = {x e R 3 | \x\ = 1, x ■ (0, 0, 1) > 0}. 

In this case, the eigenfunctions are given by those spherical harmonics which are 
odd under reflection in the {x\,X2) plane, namely, spherical harmonics 

u = Y lm = e"^P /m (cos0), \ = 1(1 + 1), 
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where —l<m<l and Z — ttt. is odd. We use spherical polar coordinates (9, <j)) 
where 

(1.4) (x±, X2, x%) = (rsin^cos^, rsin0sin</>,rcos0). 

Let us consider the case when m = I — 1. Then 

u = cie^ 1 - 1 ^ (sin 9) 1 - 1 cos 9, 
and c is a normalization factor. We have 

\H\l = c lj {(sm9) 1 - 1 cos 9) 2 sin 6 d9 = cf J ((sintf)'- 1 - (sin6>)' +1 ) 

A computation shows that 

f (sm9) l+1 d9= — l — I {sine) 1 ' 1 d9, 
Jo ' + 1 Jo 

Hence L 2 normalization requires that 

2= 1 1-21-4 2 = 1 a W i-1 \ 

' Z + 1Z-1Z-3'"! Z + l V( ? - 1 )/ 2 /' 



d(9. 



where for simplicity we took Z to be odd. This has asymptotic 
Thus, 



-2 _ , £-3/2 



2 = 4^^3/2^3/^ 

for this class of eigenfunctions. Hence there is no nontrivial lower bound of the 



form (1.1) for the hemisphere. 



In all these examples, the upper bound holds, but the lower bound fails in the 
last two examples. It is rather clear why the lower bound fails in example 3; if we 
hold m fixed and send n to infinity, then in some intuitive sense the energy in the 
eigenfunction is concentrating along lines where x is constant and these never reach 
the boundary, so most of the energy is going undetected at the boundary. More 
precisely, one can show that the semiclassical wavefront set of the family u m . n , 
with m fixed and n — > 00, is concentrating on geodesies where x is held fixed. A 
similar phenomenon is happening for the hemisphere, although here there is only 
one geodesic — the boundary curve — which is trapped, so this is a borderline 
case. Correspondingly, the lower bound is not violated as strongly in this example. 

These examples lead one to expect that the failure of the lower bound is related 
to the presence of geodesies in M which do not reach the boundary. In fact, we 
have the following result. 

Theorem 1.1. Let M be a smooth compact Riemannian manifold with boundary. 



Then the upper bound in (1.1) holds for some C independent of j. The lower 
bound in (1.1) holds provided that M can be embedded in the interior of a compact 
manifold with boundary, N, of the same dimension, such that every geodesic in M 
eventually meets the boundary of N. In particular, the lower bound holds if M is a 
subdomain of Euclidean space. 

We would like to acknowledge helpful conversations with Alan Mcintosh, Steve 
Zclditch and Johannes Sjostrand. 
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2. RELLICH-TYPE ESTIMATES 

To prove the upper bound, and the lower bound for Euclidean domains, we use 
the following Lemma which we call a Rellich-type estimate. 

Lemma 2.1. Let u be a Dirichlet eigenfunction of H . Then for any differential 
operator A, 

(2.1) / (u, [H,A]u)dg = f ^Auda. 



d 



v 



M Y 



Proof. The proof is very simple. Let A be the eigenvalue corresponding to u. We 
write [H, A] = [H — A, A] and use the fact that {H — X)u = to write the integral 
over M as 



((H - X)u, Au) - (u, (H - X)Au) dg. 
Then we apply Green's formula; since u vanishes at the boundary, only one of the 



M 



boundary terms is nonzero, giving (2.1). □ 



Let u be an _L 2 -normalized Dirichlet eigenfunction on M with eigenvalue A. To 
prove an upper bound for the L? norm of ip = d v u, we choose an operator A 



so that the right hand side of (2.1) is a positive form in ip. To do this, let us 
choose coordinates (r, y) locally near the boundary such that r is the distance to 
the boundary; this is a smooth function for r € [0, 6], where 8 is sufficiently small. 
Then we choose A — x(r)d r , where \ £ C^°(R) is identically 1 for r close to zero, 
and vanishes for r > S. The right hand side of ( |2 . 1[ ) is then precisely the square 
of the L 2 norm of ip. The left hand side may be written (after one integration by 
parts) 

(B\u, B2U) dg 



M 

where Bi are first order (vector-valued) differential operators with smooth coeffi- 
cients. This may be bounded by 

C J (Vtt, Vu) dg = CA, 

M 

where C depends on the domain, but not on A. This proves the upper bound for 
any compact Riemannian manifold with boundary. Indeed, we can prove more. 
Choosing A to be of the form Q*Qd r near the boundary, where Q is an elliptic 
differential operator of order k in the y variables, we obtain the upper bound 

(2-2) Mb"(x) ^ CXk+1 

for any integer k, and hence (by interpolation) any real k. 

Remark. Notice that this argument gives a rather explicit upper bound for 
lim sup XJ 1 \\1pjW2, namely, the reciprocal of the largest S for which r is a smooth 
function on [0, 5]. 

To prove the lower bound for Euclidean domains, M C K™, we choose A so that 



the left hand side, rather than the right hand side, of (2.1) is a positive form. We 
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choose 
(2.3) 



A = J2- 



dxi 



As is very well known in scattering theory, the commutator of this with H (which 
is minus the Euclidean Laplacian here) is [if, -A] = 2if. Hence, in this case (2.1) 
gives us 



(2.4) 



2A 



du 
dv 



Au da 



vx(^) 2 da<CU\\l 
av' 



Y 



which gives the lower bound. The equality in (2.4) was proved by Rellich |q]. 



3. Estimates for eigenfunctions near the boundary 

Let M be a smooth compact Riemannian manifold with smooth boundary. We 
choose a coordinate system (r, y) near the boundary where r is distance to the 
boundary (a smooth function for small r, say r < 6) and y is constant on geodesies 
which hit the boundary normally. In terms of these coordinates, we may write the 
metric 

g = dr 2 + hi 3 dyidyj 

and the Riemannian measure 

(3.1) dg = k 2 drdy, where fc 4 = det hij. 

Let us denote the boundary of M by Y, and write Y r for the set of points at distance 
r from the boundary, which is a submanifold for r < 5. 

Suppose that u is a Dirichlet eigenfunction for if, with eigenvalue A. It will be 
convenient to change to the function v = ku (this is equivalent to looking at the 
Laplacian acting on half-densities). Then v solves the equation 

(3.2) d 2 r v + di(h ij djv) + \ v + fv = 0, h}' J = (%) _1 
where 

/ = -k^d 2 k - k^diih^djk) 

is a smooth function on M. 

Given an eigenfunction u, we define a sort of 'energy' E(r) for each value of r. 
This is obtained formally from the energy for hyperbolic operators, with r playing 
the role of a time variable, by switching the sign of the term involving tangential 
derivatives. Let 

(3.3) E(r) = \J y (v 2 + (X + f)v 2 - h ij divd^dy. 



Lemma 3.1. For r £ [0,5], we have an estimate 
(3.4) |U(r)| < CX 



where C is independent of X. 
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Proof. From the upper bound argument in the previous section, we know that 
£7(0) = < CX. We compute the derivative of E(r): 



a 



(3.5) 9r 



— E(r) = / (d 2 vd r v + (A + f)vd r v - h^divd^ 



v 



Yr 



df 2 

+ ~-—d i vd j v+—v )dy. 
Or or ' 



Integrating by parts in the third term, and using the equation for v, we obtain 
( 3 - 6 ) ll^l^-^/ ( y2 + \ Vv \ 2 ) d y ^ C J ^(u 2 + \Vu\ 2 )k 2 dy. 

Thus, for r G [0, 8], 



E(r )=E(0)+ I —E(r)dr < CX+ I (u 2 + \Wu\ 2 )dg < CX. 
Jo dr J M 

□ 

Next, we derive an estimate on the L 2 norm of u on Y r , exploiting the fact that 
u vanishes on the boundary. 

Lemma 3.2. There exists C > 0, independent of X, such that 
(3.7) 



f s 

/ u 2 da(y) < CXr 2 for all r G [0, -] 



Remark. Given u, this bound follows, for small r, from Taylor's theorem and the 
upper bound from Section ^[ The point of this Lemma is that the estimate holds 
for r in a fixed interval, independent of X. 

Proof. Consider the L 2 norm on Y r : 

(3.8) L(r) = J u 2 k 2 dy = J v 2 dy. 

Then we have 

L '( r o) = / 2OT r dy 



and 



L"(r ) =2 (v 2 + vv rr ) dy 
Jy to 

= 2 / (yl - vdiih^v-j) - vXv - vfv) dy 
J Yro 

= 2 I (v 2 + h ij ViVj - Xv 2 ~ fv 2 ) dy 
v 2 dy -4£7(r ). 



On the other hand, from Cauchy-Schwarz we have 

(J Yro 2vv r ^ _ L'(r ) 2 



v 2 > 
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Thus we have the differential inequality for L{r): 

(3.9) L" > ^-p- - CX, 

for some constant C depending only on the manifold M. 

Intuitively, this inequality says that if L' ever becomes too big, then L" will also 
become big. This should create a feedback loop which causes L to grow rapidly, 
which we know cannot happen because we have 

s r 

2 



(3.10) / L(r)dr < u z dg=\. 

JO JM 

So we should be able to get upper bounds on V from Q and jTlcl) , which in 
turn should imply an upper bound on L. 

We use an integrating factor. After some rearrangement, (^9) becomes 

2L'L" 2(L') 3 2CXL' 

Since the left-hand side is the derivative of the quantity 

L'{r) 2 2C\ 



B(r) 



L(r) 2 L(r) 



we thus see that B is non-decreasing. 

Suppose that B(r ) and L'(tq) were positive for some < r < 5/3. Then we 
would have B(r) > for all r > r , so 

L'{rf > 2C\L(r) for all r > r . 

In particular L'(r) would be strictly positive for r > r . We rearrange this as 

(L(r) 1/2 Y = i J L'(r)L(r)- 1 / 2 > for all r > r„. 

This would give 



9 ^ 

L(r) > C'A for all r > y . 

This would contradict the bound ( |3.10| ) for large A. Thus, for A > Ao (where Ao 
depends only on M), we must have B(r) < or L'(r) < for all < r < 5/3. In 
either case 

(L(r)V2)' = ^'(r)L(r)- 1 / 2 < for all r < i, A > A . 

Since L(0) = 0, this implies (fTj]). □ 



4. The lower bound on Riemannian manifolds 

To find a lower bound on an arbitrary compact Riemannian manifold, M satisfy- 
ing the conditions of the theorem, we need to find an operator which has a positive 
commutator with H. We begin by constructing a first order pseudodifferential op- 
erator A on N which has this property to leading order, ie such that the symbol of 
i[H, A] is positive. 
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Lemma 4.1. Given any geodesic 7 in S*N, there is a first order, classical, self- 
adjoint pseudodifferential operator Q satisfying the transmission condition ( see (i) , 
section 18.2), and properly supported on N , such that the principal symbol a(i[H, Q]) 
ofi[H,Q] is nonnegative onT*M, and 

(4.1) a(i[H, Q]) > cr(H) — |£| 2 on a conic neighbourhood U 1 of 7 C\T*M. 

Proof. Recall that the principal symbol of i[H, Q] is given by the Hamilton vector 
field of H applied to q = cr(Q), the principal symbol of Q. In the case of H = 
— A g , this is |£| times the arc-length derivative along geodesies, where |£| = = 
\ note that £| is constant under geodesic flow. 

To construct q, choose a point p on 7 and a small piece of submanifold V transver- 
sal, within S*N, to 7. We next choose functions tp and b which are homogeneous 
of degree zero, so we only need to specify them on S*N. We first specify ip on V, 
requiring it to be identically 1 near p and zero outside a slightly larger neighbour- 
hood of p. We then extend t/j to S*N be requiring it to be constant along geodesies 
intersecting V, and zero on any geodesic not intersecting V. Thus ip is supported 
in a small conic neighbourhood O of 7. We define b on the support of i/j in S*N 
by prescribing it to be zero on V and to have derivative equal to 1 along geodesies, 
with respect to arc length. Then we define a, a function homogeneous of degree 
zero on T*N, by a — ijjb. Due to our geodesic condition on M, a is a smooth 
function on T*N \ 0. 

We next choose cutoff functions as follows: let x be a smooth function on N 
which is 1 on M and is supported in the interior of N, and let %2 be a smooth 
function on T*N which is 1 when |£| > 1 and when |£| < 1/2. We define symbols 
q' and q" on T*N by q' = |£|xX2a, an d q" — ~q'(x, — £). Then q" is supported 
in a small conic neighbourhood of —7 and its symbol is increasing along —7. Let 
q = q' + q"; since q is odd, we can find a classical pseudodifferential operator Q 
with symbol q satisfying the transmission condition. Replacing Q by \{Q + Q*) 
if necessary, we may assume that Q is self-adjoint. By construction, Q satisfies 



(y). □ 



We now use Lemma 4.1 to construct our operator A. For each geodesic 7 in 
S*N, we have a conic neighbourhood t/ 7 as in the Lemma. By compactness of 
S*M, a finite num ber of the U 1 cover S*M . Let A be the sum of the corresponding 
Qj. Then Lemma 4.1 implies that 



(4.2) cr(i[H,A]) > |£| 2 on T*M. 



Our next task is to get an identity for A similar to (2.1). Since A is now pseudo- 
differential and therefore non-local, it is important to take into account the fact 
that u is defined only on M. Let u denote the extension by zero to N. 

Lemma 4.2. Let u be a Dirichlet eigenfunction for H , and let A be a first order 
classical pseudodifferential operator satisfying the transmission condition. Then 

(4.3) J (u,[H,A]u)dg = 2lm J ^Auda - J (^-feda, 

M Y Y 

where c(y) = lim^oo p _1 a(0, y, p, 0). 

Remark. In ( |4.3| ), the restriction of Au to dM is taken from the interior of M; 
this is well defined since A satisfies the transmission condition and u is smooth on 
M (see §). 
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Proof. Note that u £ H 3 / 2 ~ e (N) for every e > 0, but it is not in H 3 / 2 . Conse- 
quently, we cannot unwrap the commutator [H, A] directly — since HA and AH are 
third order operators — but must use a limiting argument. Let (f>(t) be a smooth 
function which is zero for t < 1 and one for t > 2. Choose coordinates (r, y) as 
above. Notice that for < e < 5/2, (f> t (r) — (fi(r/e) may be regarded as a smooth 
function on AT, such that u e — u<f> e is a smooth function on N. 

Since [H, A] is second order, and u e — > u in H 1 (N) as e — > 0, we have 

/ (u, [77, dg = lim / (u e , [H - A, A]u e ) dg. 
Jm e ~>° Jm 

Now we can unwrap the commutator, getting 

((H - X)u e ,Au e ) - (Au e , (H - X)u e ) dg = 2 Im / {(H - X)u e , Au e ) dg. 
m Jm 

We write 

(4.4) ((H - X)u ei Au e ) = {(H - X)u e ,Au) - ((H - X)u e , A((l - 4> t )u)). 

Now we use the fact that A satisfies the transmission condition. This implies that 
Au is smooth up to the boundary of M, ie Au \ M is smooth. If we compute 
(H — X)u e , we get 

(H — X)u t = —d r ud r (j) t — ud 2 cf) e + 2uk r d r (f) e . 

The distributional limit of this as e — > is —ipSgM, and it is always supported in 
the interior of M. Therefore, 

((H - X)u e ,Au) dg — > / ipAuda 
m Jy 

where Au I M is taken as the limit from the interior of M. 



To deal with the second term in (|4.4| ), we decompose A. Let a be the principal 
symbol of A, and let y be a point on the boundary of M. Let (p, rj) be cotangent 
coordinates dual to (r, y). Then, restricted to the line rj = 0, the symbol is the same 
as that of c(y)D r for c as defined in the Lemma. Since the symbols r\i of D yi vanish 
on this line, and have independent differentials there, we can write the symbol a as 

a = c{y)p + i]ibi, 

for some symbols bi of order zero. Correspondingly, we get a decomposition 

A = cD r + B l D Vi + B', 
where Bi and B' are pseudodifferential operators of order zero. Then replacing A 



by cD r in the second term in (4.4) gives 

-/ {{H-X)u e ,cD r {u{l-(t> e )))dg -f-i / c(y)^ 2 d<j. 

Thus, to prove the Lemma it remains to show that 

/ ((H - AH, (BiD yt + B'){u{l - &))) dg -> 0. 
Jm 

Notice that, as e — > 0, the L norm of (H - X)u e is 0(e~ 1/2 ), while the L 2 norms of 
D yi (u(l — cj> e )) and of u(l — 4> e ) are 0(e 3 / 2 ). Since zeroth order pseudodifferential 
operators are bounded on L 2 , this term indeed goes to zero. This completes the 
proof of the Lemma. □ 



10 



ANDREW HASSELL AND TERENCE TAO 



Now we use (PI) and Q to prove the lower bound. First, we deal with the 
left hand side of (4.3). Using (|4.2|), we may write 



i[H,A] = H + B*B + R-L+ R 2 , 

where B and R\ are first order pseudodifferential operators, and R2 is second order 
with kernel smooth on M x M. Thus, we can estimate from below 



(4.5) 



(% i[H, A]u) dg > A + / (u, (Ri + R 2 )u) dg 



M 



> \-C\\u\\ Hl \\u\\ L 2 =X-CX 1/2 . 



We want to show that this is no bigger than CH^Hl- To do this, it is sufficient to 
show that 



(4.6) 



\\Au lY \\l HY) <CX + C\\ 



since by Cauchy-Schwarz we can then estimate the right hand side of (4.3) by 



(ip, Au) da < - 1 



e\\Au\ 



Y\\L 2 (Y) 



use (L6) and absorb the ||^4w|| 2 term in th e lef t hand side. This, to finish the proof 
of the lower bound, it is enough to prove (4.6). 



To do this, we require the following Lemma, which is due to Boutet de Monvel 
(lj and Visik and Eskin || . 

Lemma 4.3. Let B be a classical pseudodifferential operator on N of integral or- 
der k > — 1 satisfying the transmission condition. Let B r denote the operator 
from L 2 (Y r ) to L 2 (Y) given by B r (f r ) = B{f r ® 8y r ) \Y. (For r = 0, we take 
lim e ^o(B(/o ® <5y))(e, •); this is well defined for operators satisfying the transmis- 
sion condition.) Then we have an estimate 

(4.7) \\B r \\ LHYr )^LHY) < Cr- k -\ r e (0,6]. 

For k — — 1 we have B r uniformly L 2 bounded down to and including r = 0; indeed, 
in this case, if we identify Y r and Y , then {B r } is a uniformly bounded family of 
zeroth order pseudodifferential operators on Y , for r G [0, S\. 

Proof. The proof of this is at least implicitly contained in [Q, in the discussion 
before Theorem 18.2.17, but we give a self-contained proof here which is based on 
that in Q. First, since ^ k (N) is invariant under adjoints, it is sufficient to bound 
B mapping from L 2 (Y) to L 2 (Y r ), which by an abuse of notation we shall denote 
B r for the duration of this proof. 

Let the kernel of B be given as an oscillatory integral near r = r' = by 

B(r,y,r',y') = J e^^e^^^b^y,^')^!^'. 

The kernel of B e is obtained by restricting this to r = e, r' = 0. Then we wish to 
obtain a uniform L 2 bound on r k+1 B r . Multiplying the kernel by r k+1 is the same 
as applying D 1 ^ 1 to the phase. Integrating these derivatives by parts reduces the 
order of the symbol to —1. Thus, it is enough to consider the case k = —1. 

Now consider the operator B r acting on a fixed function v G L 2 (Y) . Equivalently, 
we can look at B acting on vSy- This is the distributional limit of B acting on 
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ve 1 (f)(r/e), where <f) £ C^°(M) has integral one and is supported in \r\ < 1, say. 
The result is 

(4.8) J e^e^'bir, y, 0v{?)$(<&) <%i <• 
Thus, the symbol of B r is 

(4.9) lim / e^bfay^Z'Me&dS!. 



If b in (4.9) were a symbol of order —2, then the limit in fl4.9|) is 



(4.10) e^b(r,y,^,ad^. 



which is a symbol of order —1 in the £' variables, depending uniformly on r (ie, 
with uniform estimates on all seminorms). For b of order —1, we choose a cutoff 
function x(£i)j where x is supported in |£i| < 1 and is identically one near £i = 0, 
and write b = x(£i)^ + b' . Then x(£i)& is a symbol of order —1 in the £' variables, 
uniformly in r and £i (since it is compactly supported in £i), so this integral gives 
a symbol of order —1 uniformly in r. Finally, b' is supported away from £j = 0, so 
we may write using Taylor's theorem 

n 

&' = (i-x(a))( c (r,y)er 1 + E^)' 

where 6j are symbols of order —2. Using the statement above about symbols of 
order —2, we see that the ^bj terms give symbols of order 0, uniformly in r. On 



the other hand, the £f 1 term can be calculated explicitly. In the integral (4.9), 
with b replaced by (1 — x(£i)) c ( r ; y)£i j the integral can be moved into the upper 
half complex £i plane provided e < r, and then the result is 

c{r,y)[J e^(l-x(Ci))rti + / e^f 1 ^), 

|£i|=l, 0<argfi<7r 

which is a multiplication operator depending smoothly on r. This completes the 
proof. □ 

Let H^ 1 denote the inverse of the Laplacian, with Dirichlet boundary conditions, 
on N. Then, in a neighbourhood of M, H^ 1 is a classical pseudodifferential operator 
of order —2 satisfying the transmission condition. Since Hu = Xu + ipdy , we find 
u = H- 1 (\u + i/)6 Y ). Thus, 

Au = AH- x (\u + ip5 Y )- 



By Lemma <L3, AH~ 1 (i/j8y) \ Y is bounded by C||-0||2, which satisfies ( |4.6| ). Let us 
consider \AH~ l u. To do this we break up u into a 'close' and 'far' part. Since 
the length scale is A -1 / 2 , let us use the functions (j) e and u € as above, b ut c hoosing 
the value e = A -1 / 2 . Consider first AAiJ _1 (lt(l — (f> e ))- By Lemma |4.3| , we get 
an estimate independent of r on the operator norm of (A£T -1 ) r , and u(l — <fi e ) is 



supported in r < 2A 1 / 2 , so we can estimate using (3.7) 



\\XAH- 1 (u(l - e ))||a < C-^ / rA 1/2 dr < CA 1 / 2 , 
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which satisfies (4.6). Finally, consider the term XAH 1 u IL . Writing u e — (A 2 H 2 u) 
we obtain 

AAff- V = \~ l AHUz - X^AH^lH 2 , <f> e ]u. 



The first term is estimated via Lemma 4.3 and ( p.7| ), and using smoothness of the 
kernel of AH at (r, y, r' , y') for r = and r' > S/3, by 

A -1 C f S/3 r-*\V 2 rdr + 0(lj) ^CfX^ + X- 1 ) 

which satisfies (4.6) for large A. The terms coming from the commutator are treated 
similarly. 



This completes the proof of Theorem 1.1 



5. AN ALTERNATIVE PROOF OF THE LOWER BOUND 

There is an alternative way to prove the lower bound, which depends on a trick 
due to Morawetz, Ralston and Strauss ||. First, we recall some facts about spheri- 
cal harmonics. Let Aj»-i denote the Laplacian on the (n— l)-sphere. Then — Agn-i 
has eigenvalues k(n + k — 2), k = 0, 1, 2, . . . . Let the corresponding eigenspace be 
denoted elements of 14 are called spherical harmonics. We recall that for ev- 
ery (j> G Vfc, the function r k cf> (thought of as a function on K™ written in polar 
coordinates) is a homogeneous polynomial, of degree k, on R™. 

The alternative proof goes as follows. We use the symbol a of the operator A 
constructed in the previous section, restrict it to the cosphere bundle of N, and 
expand in spherical harmonics. Note that since a is odd, and spherical harmonics 
are even or odd with k, all the terms with k even vanish. Thus we have 

oo 

a\S*N = J2 fcmfo jfi)> M*, ') e V k (S*N). 
i=o 1^1 

Since a is smooth, this sum converges in C°°, and certainly in C 1 . Therefore, there 
is a K such that the operator A' with symbol 

a ' = ^2i+i(a;, j7|) 
i=o 

also has positive commutator with H. However, following j^], we can turn A' into 
a differential operator P of order 2K — 1, by letting 

R7-1 



i=0 

The symbol p is a polynomial on each fibre of T*N, so we may take P to be a 
differential operator. Moreover, the symbol of i[H,P) satisfies 

a(i[H,P]) = \£\ 2K (<7(i[H,P])\ lil=1 ) > c\£\ 2K for some c> 0. 

We now apply apply the Garding inequality to Q = i[H,P]. We need to apply 
the form which is valid for u £ H K (M) as opposed to Hq (M) which is the standard 
situation. This gives 

. K-l 

(5.1) / {u,Qu)dg > c\\u\\ 2 h k (m) - c(\\u\\ 2 L 2 (M) + ll^u||^K-i/ 2 - fc(r) ), 
Jm k=o 
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where c is a positive constant depending on P and [M, g) . This can be derived from 
the usual Garding inequality by writing u = v + w, where v solves the Dirichlet 
problem (P + a)v = with the same Dirichlet data as u. (Here a is a sufficiently 
large constant.) We then get an estimate on the H K norm of v in terms of the 



boundary norms of u appearing in (5.1). On the other hand, we can apply the 
standard Garding inequality to w. 
Thus, ( pH] ) and ([Tl]) combined yield 

K-l 

2 i f~i \ \\ fik 



\ U \\H K {M) 



<C + CJ2 ||^u||^- 1/2 - fe(y) + I / (i>,Pu)da 

k=Q ^ Y 



We can estimate the left hand side from below by a constant times X K . Consider 



the right hand side. Let us write u — k~ 1 v, where k is as in (3.1), so that v satisfies 
([3~2]). Then Pu = Pv, where P — P o k is a differential operator of order 2K — 1. 
Since k is smooth, we obtain 

A'-l „ 

(5.2) \ K <C + CY, \\drv\\ 2 H K-i/2- HY) +C\ / (iP,Pv)da\. 

k=0 ^ Y 

Since v = at Y, and we are interested in Pv \ Y , we may assume that P = P'd ri 



where P' has order 2K — 2. Using (3.2), we may replace d^v by —(A — f)v 
d Vi (h lJ d Vj v) repeatedly, until only d r dyV terms remain. Thus we have 



K-l 

3=0 



where Pj is a differential operator on Y of order 2(K — 1 — j), independent of A. 
Hence (5.2) becomes 

K-l K-l 

(5.3) X K <C + CY / X^WiifxK-i/z-H^+CY,*] / (^,P^)da\. 



k=0 j=0 



Using the upper bound estimate {2/2) on the sum over k and for all terms in the 
sum over j with j < K — 1 , we find 

(5.4) \ K < c + c(i + x K - 1/2 ) + cx^uwl* + cx K - 1/2 u\\ L 2. 

Finally, we estimate the last term in (5^) by eX K + e~ 1 C 2 X K ~ 1 HV'lll) an d absorb 



the eX K term in the left hand side. This proves the lower bound. 

6. Further examples 

Here we briefly discuss further examples which illustrate how the presence of 
trapped geodesies in the interior of M affects the behaviour of eigenfunctions at 
the boundary. 

Example 5 — the spherical cylinder. Consider the manifold M C S 2 given in 



terms of coordinates {9, 4>) as in (1.4), by 



M = {(6,<P)\^-a<0<^+a, q^eS 1 }, < a < |. 

Let y = — it/2. Here, the equator at y = is a trapped geodesic, as is every nearby 
geodesic, so this is an extreme case where there is a nonempty open set of periodic 
geodesies. In this case, there is a sequence of normalized Dirichlet eigenfunctions 
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ui where the lower bound fails spectacularly, namely there is exponential decay of 
(6.1) \ml<Ce~ d 

This follows by separating variables, u = e ^v(y), and reducing to a one dimen- 
sional problem where l~ l plays the role of a semiclassical parameter. The positive 
curvature of the sphere manifests itself as a potential function with a nondegen- 
erate minimum at y = 0. Well known estimates (see e.g. chapter 3) show 
that there are normalized eigenfunctions it; which localize near y — as I — > oo. 
More precisely, for any 8 > and differential operator P there are constants c > 
and C < oo such that \Pui\ < Ce~ cl in the set {|y| > S}. If we then apply our 
upper bound argument from Section wh ich involves only the gradient of the 



eigenfunction near the boundary, we obtain (6.1) 



Example 6 — a finite hyperbolic cylinder. Let M be the Riemannian manifold 

M = [-a, a] y x Si, g = dy 2 + (cosh y) 2 d(j) 2 . 

Again, there is a periodic geodesic at y = 0. However, this is the opposite extreme 
to the example above. Now M has negative curvature, the geodesic is unstable, and 
it is in fact the unique trapped geodesic on M. (Other geodesies may be trapped 
as t -> ooori^ — oo, but not both.) Using results of Colin de Verdiere and 
Parisse one can show that there is a sequence of normalized eigenfunctions it;, 
I = 1,2,3, ... , such that 

(6-2) r^r<\Ml< CX 



log A log A 

Thus, the lower bound is violated, but just barely. Both this example and the pre- 
vious one are based on sequences of eigenfunctions which microlocally concentrate 
at a periodic geodesic. Since the geodesic is in the interior of the manifold, the 
eigenfunctions are concentrated in the interior, and so are not as large near the 
boundary as one would normally expect. 

Example 7 — Neumann eigenfunctions. One might hope to prove bounds of the 
form 

c < \\u iy \\l 2 (y) < C 

for normalized Neumann eigenfunctions u. The following examples show that nei- 
ther the upper nor lower bound is valid in general. 

Consider the unit disc, as in Example 1. There is a basis of Neumann eigenfunc- 
tions of the form 

where Ck, n is a normalization constant and j' k n denotes the kth zero of the derivative 
of the nth Bessel function J n . The eigenvalue for this eigenfunction is A = (j' k n ) 2 . 
Denoting the boundary value of this eigenfunction by Xk,m Rellich's identity (§□]) 
gives us 

(6.3) 2A=-/ (u,d 2 u)d9= [ (u,(\ + d 2 )u)d9=(\~n 2 )\\xk,n\\i 

Js 1 Js 1 
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This shows that there is a lower bound for ||xfe,n||2- However, the first positive 
critical point of J„, j[ n , is known to have asymptotic 

fin ~ n + en 1 / 3 , n — > oo. 



So fixing fc = 1 and sending n — > oo in (S.3) shows that there is no upper bound 
for ||Xfe,n||2- 

We remark that it is not hard to derive a lower bound for convex Euclidean 
domains from Rcllich's identity. On the other hand, the argument of Example 5 
applies equally to Neumann eigenfunctions, so the lower bound is not generally 
valid. It is not clear to the authors whether the lower bound holds under the 



nontrapping assumption of Theorem 1.1 



Example 8 — vector fields are not enough. One might wonder whether, on 
an arbitrary compact Riemannian manifold, with no trapped geodesies (that is, 



satisfying the conditions of Theorem 1.1), one could choose a first order differential 
operator A whose commutator with H had a positive symbol. In this section, we 
show via an example that this is impossible in general. First we analyze what it 
means for a vector field to have a positive commutator with H . Let the symbol of A 
be ai{x)^i. The Hamilton vector field of H is £,id Xi , so having a positive commutator 
requires that 

(6.4) ^.|g >0 for |£|^0. 

Thus, the matrix d Xi aj must be positive definite. 




Now consider the 2-manifold M with boundary illustrated in the figure (with 
corners smoothed so that it has smooth boundary) — in the figure, the top and 
bottom dashed lines, and the leftmost and rightmost dashed lines, are identified. 
It is clear that M satisfies the geodesic condition of the Theorem, so the lower 
bound is valid for it. Assume that there is a vector field A on M having positive 
commutator with H. Notice that the two points p and q are such that there are 
three geodesies from p to q: one in the direction e\ + e2, one in the direction —e\ 
and one in the direction — e2. Write A = a\e\ +a2C2- Then (6.4) implies that a\ is 
increasing in direction e%, and a2 is increasing in direction e^- Hence a\(p) > a\(q), 
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and a,2(p) > »2(g/)- On the other hand, a\ + a 2 is increasing in direction e\ + e%, so 
this yields a\(p) + a 2 {p) < a\{q) + 012(9), which is a contradiction. Hence no such 
vector field A exists. 
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